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Abstract 
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We study branes residing in infinite volume space and of finite extent in the transverse directions. We 
calculate the graviton propagator in the harmonic gauge both inside and outside the brane and discuss 
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J> , its dependence on the thickness of the brane. Our treatment includes the full tensor structure of the 
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CO ' propagator. We obtain two infinite towers of massive modes and tachyonic ghosts. In the thin-brane 

o : 

limit, we recover four-dimensional Einstein gravity. We compare our results to similar recent results by 
^ . Dubovsky and Rubakov. 
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I. INTRODUCTION 



The weakness of the gravitational force has been successfully explained by postulating the 
existence of extra dimensions 



DDDQ . The effect of the extra dimensions is a high-energy 
modification of Newton's Law of gravity due to the tower of Kaluza-Klein modes. When the 
extra dimensions are of infinite volume, light Kaluza-Klein modes may dominate even at low 
energies Ifil, LZl 1^11 - Thus, unlike with finite- volume extra space, Newton's Law is modified at 
—lily L ee distances flUHQHUHy Q. Dvali a nd Gabadadze Q showed 
that this is not the case if the infinite space in which the brane lives has dimension D > 5. 
They studied a three-brane of the 5-function type and showed that the graviton propagator 
has a four- dimensional momentum dependence on the brane even at low energies. This feature 
is expected to persist if the brane is of finite thickness ("fat") in the transverse directions for 
phenomenologically relevant values of the momentum. For extremely low energies, a fat brane 
should lead to a higher-dimensional behavior of the propagator. This was discussed qualitatively 
in Q. 

Here, we present a quantitative study of a fat brane in infinite volume extra space. By 
linearizing gravity in the harmonic gauge, we arrive at an explicit expression for the graviton 
propagator. First, we obtain the propagator for the trace of the metric field over the transverse 
directions. The trace is a scalar field from the four-dimensional brane point of view. This scalar 
then contributes to the four-dimensional graviton propagator as a source, in addition to the 
matter fields. This complicates the tensor structure of the graviton propagator which becomes 
momentum dependent. We explicitly obtain the solution for the graviton propagator and analyze 
its momentum dependence and pole structure. We find two infinite towers of massive modes and 
tachyonic ghosts. In the thin-brane limit we recover four- dimensional Einstein gravity on the 
brane. 

Our discussion is organized as follows. In section [HI a brane-bulk action is considered which is 
similar to the action of ref. |7[ but generalized to allow for a brane of finite thickness in the bulk. 
In section IIII1 we solve the linearized Einstein field equations in the harmonic gauge and obtain 
an explicit expression for the graviton propagator. In section IIV[ we analyze the pole structure 
of the graviton propagator and compare our results to the Dubovsky-Rubakov model ^6]. In 
section we analyze the momentum dependence of the graviton propagator both on the brane 
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and in the bulk. We finally conclude with a summary of our results in I VII 



II. ^-DIMENSIONAL FAT BRANE MODEL 

We are interested in the dynamics of a 3-brane in a D-dimensional infinite space. The 3-brane 
is allowed to have finite thickness in the bulk-space with extent governed by the density function 
(T\{y). The action is similar to the one discussed in ref. [3], 

S = M D - 2 J d 4 xd D ~ 4 yV=9 ^ {D) + M 2 J d 4 xd D ~ 4 y^ ^{v)^ + S matteT (1) 
where qab is the .D- dimensional metric which generates the D-dimensional Ricci scalar 1ZS D \ 

— (4) _ 

whereas 1Z is generated by the four- dimensional metric g which is the induced metric on the 
slice y = const.. Capital Latin indices run over D-dimensional space-time (A, B = 0, 1, 2, D), 
Greek indices run over the four-dimensional brane worldvolume spanned by coordinates x^ (/i = 
0, 1, 2, 3) and lowercase Latin indices run over the extra space spanned by y m (m = 4, 5, D). 
M is the .D-dimensional Plank mass. S matter is the (unspecified) matter action giving rise to the 



fat brane configuration. We set y = \y\ = J y\ + y\ + ... + y%_ 4 and a^{y) is a smooth function 
of width 1/A approximating a 5-function. The mass scale M is related to the four- dimensional 
Newton constant. In general, M will depend on M, but here they will be treated as independent 
scales. 

For explicit calculations, we will choose a step-function form of the density <j\, 

(D-A)A D - 4 

= { - 1 0(1/A - y) (2) 

where u n is the surface area of the unit n-dimensional sphere. The careful reader may wish to 
smoothen the step-function first and then take the limit in which <ta becomes discontinuous. Our 
results are not altered. In the limit A — > oo, the density o"a approaches a 5-function: <J\(y) — > 
5 D ~ A (y). This continuous distribution of the 3-brane may be thought of as the continuous limit 
of a discrete set of the four- dimensional hypersurfaces (infinitely thin 3-branes) discussed in Q|. 
The Einstein field equations are 

M D - 2 Gf^,y m ) + M 2 a A (y)Gf B (x^y m ) = T AB (x^y m ) , (3) 
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where G^ is the .D-dimensional Einstein tensor and G AB only has brane worldvolume compo- 
nents. Expanding around a flat background, 

9ab = Vab + h AB (4) 

the first-order Einstein equations are as follows. The transverse components give 

2d A d n h An - d n d n h A - d A d A h n n = (D - A)(d c d D h CD - d c d c h%) (5) 

The mixed components give 

d A d A h an = d a d A h An + d n d A h Aa - d a d n h A (6) 

and the brane worldvolume components imply 

M D ~\d a d A h pA + d p d A K A - d A d A h a p - d a d p h A - Vaf) {d A d B h AB - d B d B h A ) ) 
+M 2 a A (y) {d a d v hp v + d p d v h au - dvVhpc - d a dph v v - VaP (d»d» V - d^K) ) 

= T aP (x»,y) (7) 

where we have chosen a matter source described by the stress-energy tensor T^ u whose transverse 
components vanish (T mn = 7j, n = 0). Indices are raised and lowered by the flat metric tensor 
Vab- 

To solve the field equations, we shall choose the harmonic gauge, 

d A h AB = U B h A . (8) 

We obtain from eqs. © and (jEJ), respectively, 

(6 - D)d A d A h n n = (D-A)d A d A h» (9) 
d A d A h ma = (10) 



so we may set 



h ma = (11) 
(D - 6)K + (D - 4)h* = (12) 
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Then the brane worldvolume components of the Einstein equations can be written in the following 
form: 

- M D ~ 2 d A d A (Kp - \<n a(3 h B B ) + MVa {-dyVhap + d a dph n n - fraf, d»dM - K)) 

= T aP (x^y) (13) 

Performing a Fourier transform in the brane worldvolume coordinates x M and multiplying by 
an arbitrary conserved stress-energy tensor T' al3 , which for simplicity is assumed to have no 
^-dependence, we obtain 

(m d ~ V + dnd n ) + MV A p 2 ) Kpf'^ = fapW, y)f"*P 

+ \f'» (mV a p 2 (~K ~ K) + M D ~ 2 (p 2 + d n d n )l4f) 

where the Fourier transformed, D-dimensional d'Alembertian is 8a9 a = —d n d n — p 2 with p 2 = 
Pq — p 2 the worldvolume Minkowski four-momentum. In the next section, we shall solve this 
equation for the graviton propagator. 



III. GRAVITON PROPAGATOR 



In general, the spread functions of the brane and the matter source are different. However, 



it was argued by Dvali, et al. |25( that the two spreads coincide at lowest order with correction 
terms suppressed by factors o(M/M). We shall therefore adopt a source stress-energy tensor of 
the form 

T aP (x^y)=T al3 (x^a A (y) (15) 

in the explicit calculation of the tensor structure and momentum dependence of the graviton 
propagator. Taking the trace of eq. (jlHj) . we obtain 

- M D - 2 (p 2 + d s d°)K + ^| MV A p 2 K = f: a A (y) (16) 

where we used eq. ()12|) to express in terms of h™. This is an equation for the field (trace 
over transverse directions of the metric field), which is a scalar from a four- dimensional point of 
view. The solution is obtained on the brane and in the bulk in terms of the Green function to 
the wave equation, 

(m d - 2 (p 2 + d s d s ) - (A - 1) M 2 pV A ) Q x (p, y) = a A (y) (17) 



as 



(D-2) a ' ^ ' P-2) 
After some algebra, we obtain a spherically symmetric solution expressed in terms of Bessel 
functions as 

1 / 1/1 \ (D ~ 6)/2 \ 

v) = ~ mV(a - 1) I 1 " A U J ^/.(p/^eiMM.^) j ( 19 > 

inside the brane (y < 1/A), and 

i a a / 1 \ (D - 6,/2 



5a(? ' ^ = ' mV(A-D t J W(WA)B^W (20) 

in the bulk (y > A), where 

B x = k x I (D ^ )/2 (k x p/A)H$_ 6)/2 (p/A) + I {D _ 6)/2 (k x p/A)H^_ 4)/2 (p/A) (21) 
and we have introduced the constant k x given by 

k>-(\ I) {D ~ 4)AD " 4 ^ l~fA 1) P~4)A g - 4 M 2 

^ " ( A ~ 1} u D _ 4 - 1 ~ (A - 1) Wd _ 4 mD _ 2 (22) 

To obtain the graviton propagator, we will also need the Green function which is the solution to 
eq. (JT7j) when A = 0. Notice that when A = 0, eq. (fT7|) turns into the wave equation for a scalar 
field in the thin-brane limit. It is derivable from a scalar field action. Explicitly, 

if 1/1 \ (D ~ 6)/2 \ 
Go( P ,y) = I 1 + -p i^j H$_ m (p/A)J {D _ 6)/2 (Kpy)\ (23) 

inside the brane (y < 1/A), and 

(D-6)/2 

^ (-) 

P 

in the bulk (y > 1/A), where 

S = KH$_ 6)/2 (p/A)J iD „ 4)/2 (K P /A) - HU_ 4)/2 (p/A)J {D _ 6)/2 (Kp/A) (25) 

and (see eq. 



1 / 1 A ( D ~ 6 )/ 2 
fi)(P.y) = -^g (^J J (D - 4 )/2(^/A)i/ ( ( l 6)/2 (ra) (24) 



(D-4)A D ~ A M 2 {D-A)A D ^M 2 



= ~ k '' = 1 + „ C J " U J, M°-> (26) 
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Notice that inside the brane, Go{p,y) oscillates rapidly over the transverse width of the brane. 

We are now ready to deduce the full graviton propagator. To this end, let us massage eq. (fTljl 
into the form 



M D ~ 2 (p 2 + d n d n ) + M 2 p 2 a A (y) ) { hafiMf** + 3p _4) Kfa v)Tl v 

r^f^-if/^U A (y) (27) 



The solution for the graviton propagator is readily obtained in terms of the scalar propagators, 

Mp, y )f«# = jr^r ** - ir^} g ( P , y ) + x D Z% f ffi s &> ^ ( 28 ^ 

where we used eqs. (fTTj) and (JTHJ). 



IV. POLES OF THE GRAVITON PROPAGATOR 

Next, we analyze the pole structure of the graviton propagator. We then compare the results 
of our model with that of Dubovsky and Rubakov j^ . 

A. Our model 

Using the expressions for Go(p,y) and (fE?|) for Gx(p,y), the graviton propagator (J28j) 
inside the brane (y < 1/A) can be written in the form 



r i / i \ (o-6)/2 



ATp 2 



\fafif** ~ \f:ff} ^J {D - m {Kpy) + \f:ff^I {D _ m {k x py) 



(29) 



For convenience, we have separated the term that corresponds to the tensor structure and mo- 
mentum dependence of the four-dimensional graviton propagator. To study the pole structure, 
we shall introduce the average value of the graviton propagator over the transverse directions 
of the brane (see for problems associated with the definition of observables on the brane) 
defined by 

hlT^P) = [ d D - 4 ya A ( y )h a p(p, y ) (30) 



Integrating (J29j). we obtain 



ft=r (p)t- = {f„ s r- - ^} + - if WKp2/A)2|i _ ^ ((tp/A)| 

_ ifaf/p (-P - 6) /o-j\ 

6 a * (MA: A pVA) 2 [l + ^(WA)] 
from which we may easily deduce the pole structure of the graviton propagator inside the brane. 
The above expression is valid for D > 6 (for D = 6, we obtain logarithmic corrections, but the 
results are similar and will not be explicitly discussed here). The functions that appear in the 
denominators in (j31j) are 

Z J{D-4)/2{Z) Z I(D-4)/2{z) 

for D > 6. The poles of the propagator are solutions to the equations 

/xo(«p/A) = l , /j,\(k\p/ A) = — 1 (33) 
Using (122) and the Bessel function identity 

zJ u -i(z) + z.J v+l (z) = 2vJ v (z) (34) 

for v = (D — 6)/2, it is easily shown that the solutions to fio(z) = 1 are the roots of J I/ _ 1 = 
J(D-8)/2- As is well-known, there are infinitely many zeros for v > 0, i.e., D > 6, which is the 
case we are considering here. We shall denote them by Zj, 

J(D-s)/2(zj) = , J = 1,2,... (35) 



We therefore obtain an infinite tower of massive poles with masses given by 

A 2 

K 

Similarly, the condition (jL\(z) = —1, together with the Bessel function identity 



m) = z) -J (36) 



z/„_i(z) - zl u+l {z) = 2vl v (z) (37) 

and the relation I v {z) = e^ nu ^ 2 J u (iz), lead to a tower of tachyonic poles with masses given by 

A 2 

m % = ~ z ) ^2 (38) 
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To obtain the behavior of the propagator near a massive pole, observe that 



J(D-8)/2(z) J{ 



i - «,w = = - 7 D - 8)/ f j ; ( , - + <<(* - ,,) 2 ) (39) 

Using the Bessel function identity 

zJ' v _ x {z) = {v- 1) J^(z) - zJ v {z) (40) 

together with (jSljl. we deduce 

1 -M^) = ^ L ^(^ 2 -^) + ... (41) 

near z = Zj. It follows that the graviton propagator on the brane (j31j) behaves as 

~ h ^ r{p)fla p „ ^f apfla p _ 2^)~ 6) V*? (42) 

near the massive pole p 2 = m 2 . Similarly, near the tachyonic pole p 2 = mL, we obtain 

1 ~ -m 2(D-6) 2 /z^ 
h^ e (p)T' al3 ~ --T^T'f ^ } -Ul (43) 

The minus sign of the residue of the tachyon implies that the tachyon is a ghost. 

Notice that both the massive modes (jS^jl and the tachyons (|S5|l are expressed in terms of the 
same mass scale parameter p c , where 

, A 2 A 2 M D ~ 2 , . 

Pc 2 ~ 72" ~ =2 44 

In the thin-brane limit (A — > oo), we have p c — > and the infinite twoers of massive modes and 
tachyons turns into continuous spectra. The form of the propagator in this limit is easily deduced 
from eq. (j31)l . For momenta away from the critical scale (\p\ ^> p c ), the two terms in (|3*T|) that 
give rise to the massive and tachyonic poles become vanishingly small and we are left with 



hir e (p)T iaP ~ {tap** - \f:ff) * (45) 

l 2 ) M p 2 



recovering four-dimensional Einstein gravity. 
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B. The Dubovsky-Rubakov model 

It is interesting to note that similar results have been obtained by Dubovsky and Rubakov 
using a slightly different model. In order to directly compare our results with theirs, we shall 

n 

assume that the spread function (denoted by f 2 (y) in [26]) is given by eq. (J2J). Then the Einstein 
field equations proposed in |26J] can be written as 

F{nW)G%{x>, y) + M 2 a A (y) j d D -*y>a A (y') Gf B (x», jf) = T AB {x^ y) (46) 

to be compared with the Einstein eq. (jBJ) in our model. In eq. ( |46jl . the four- dimensional Einstein 
tensor only has brane worldvolume components (i.e., G^ B = 0) and the form-factor T ~ M D ~~ 2 
at low energies. Also, the matter source on the brane will be assumed to have only space-time 
components and a spread function same as that of the brane, 

{x)cr A (y) (47) 

where T M ^(x) is conserved in the four- dimensional sense (cf. eq . (J15|) in our model). The inverse 



width A of the spread function is assumed to be A ~ M in 



2611 to be contrasted with our model 



in which A ~ M, since it coincides with the inverse width of the brane Q|. 

Working as in section |Hj we linearize the Einstein equations and obtain the graviton propa- 
gator in the form 

Mp, y)f'^ = I \f, v f'^ - Iffi^ Q 1 (p, y) - ^Lffi^ip, y) (48) 

where we multiplied by the arbitrary stress-energy tensor T' to absorb the longitudinal part 
which is not gauge-invariant. It is given in terms of the Green function which satisfies eq. (JT7j) 
for A = 1, 

m d ~ 2 (p 2 + d s d s )Gi( P , y) = My) (49) 

i — j 

(denoted by Df in |26j). The denominators are 

c = \ + M 2 p 2 g 1 , c* = i-M 2 P 2 g 1 (so) 

where fijis the average of Q\ over the spread function (defined as in eq. (fHUj) and denoted by 
Df f in |2£|). Explicitly, 



01(P) = -=S^/(P/A) , f(z) = \ 
MA 2 z 2 



171 



l__L (jD _4)if ( «_ 4)/2 (z)J (D _ 4)/2 (z) 



(51) 
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«> - (j - 4 »rr (52) 



where we introduced the function f(z) for convenience and the scale k, which coincides with our 
earlier definition (ffilj) in the large A limit, 

The poles of the propagator (}4*%j) are the zeros of C and C*. They can easily be seen to correspond 
to small z, therefore we may approximate C ~ 1 — k 2 f (0)p 2 / A 2 , whose root is 

A 2 M D - 2 

m 2 ps - - - ~ — 5 , (53) 

which is a massive pole. Similarly, the root of C* is a tachyonic pole 

2 2 , x 

m t ~ — m ~ o (54) 



Notice that the mass scale is similar to the mass scale of the poles in our model (|44|) . although in 
this model only one pair of poles is obtained instead of the infinite tower we found in our model. 
This scale matches the one found in 2(| if we set A ~ M, in which case m ~ M 2 /M. 



V. MOMENTUM DEPENDENCE OF THE GRAVITON PROPAGATOR 

Having understood the large A limit, we now turn to a study of the momentum dependence 
of the graviton propagator keeping A finite. By introducing the width 1/A, we have added a 
scale to the theory in addition to the mass scales M and M. It follows from the explicit form 
of the propagator that the relevant scales are A and A/k, where k ~ k\ ~ k is a dimensionless 
parameter given by f!26[) or 122(1 . Phenomenologically, one expects A ~ M and M M. So 
we shall restrict attention to momenta that are well below the scale A (p <C A). This range is 
divided by the scale given by eq. (|4*4"|) into a small momentum (p p c ) and a large momentum 
(p 3> Pc) regime. Qualitatively, one expects four- dimensional behavior of the graviton propagator 
for large momenta and D-dimensional behavior for small momenta [7|. We wish to study this 
behavior quantitatively. 

For small momentum, p p c , we have 

Go(p,y)*Gx(p,y) (55) 



n 



as can easily be verified from eqs. (|2Uj) and (J24j) in the bulk and eqs. (119)) and ()23|) on the brane. 
The resulting tensor structure of the graviton propagator (}2~%j) is 



K P { P] y)P# « jf^f** - ^l^f^f/} g (p, y) (56) 



In the bulk, we deduce from 

, -. x (X>-6)/2 

&(kv)~*(-J ^ ( ( l 6)/2 (py) (57) 

which is the propagator for a D-dimensional scalar field. Therefore, the graviton behaves as a 
.D-dimensional field in both its momentum dependence and its tensor structure in the bulk. 
On the brane, after averaging over its transverse width, eq. (joljj) yields in the regime p<j) c 



M p l I \ u Z J 



where we used eqs. fl2H)) and (J3"U|) . It is easy to see that the 1/p 2 pole vanishes. The first 
non-analytic term can be found from the expansion for small argument 

H M = -|r M (|) -(! + ...) + (£)"■■ (§) + ... (59) 

for integer z/, where the dots represent higher-order and analytic terms. Applying this to eq. ()58|) . 
we obtain 

h*p»(p)T*# ~ jf^f** - ^i—ft^? J (|) ln(p/A) (60) 

exhibiting D-dimensional behavior. Similar conclusions may be drawn for the trace h™ in the 
small momentum regime (p <C p c ). 

In the large momentum regime (A ^> p ^> p c ), the results are similar to those in the large A 
limit, which we discussed in the previous section. In this regime, the scalar Green functions are 
related by 

Go(p,y)K{\-l)Qx(p,y) (61) 

to be contrasted with the relation (|55|) in the regime p <C p c . Thus the tensor structure of the 
graviton propagator (J25)l becomes 

h aP (p; y)f'^ « {f aP f^ - \TZf£) g (p, y) (62) 
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exhibiting four- dimensional behavior. Inside the brane, we deduce from (J23j) 

~ h B r{p)fla , = Ufn* _ IfafA 1 + Q ((pc/p)2) ( g 3) 

I A ) M p 2 

exhibiting the distance dependence of Newtonian gravity with the tensor structure of four- 
dimensional Einstein gravity. This is in agreement with our earlier conclusion ()45|) in the large 
A limit. In the bulk (y > 1/A), we deduce from ()24|) 

i / 2 \ (-D-6)/2 ( 1 1 / 1 \ ( £, ~ 6 )/ 2 

MP, y)T- - =^ Qk) { V« - i^} (1) <_ 6)/2 W 

(64) 

Therefore, the propagator vanishes in the thin brane limit (A — > oo). These results for the 
momentum dependence of the tensor structure of the graviton propagator are in agreement with 
the qualitative suggestions of Q. 



VI. CONCLUSION 



The main objective of the present work was to analyze quantitatively the gravitational effects 
on a brane of finite extent in the transverse directions ("fat"). We obtained an equation for 
the graviton propagator which we then proceeded to solve in two steps. First, we obtained the 
propagator for the trace /i™ over the transverse directions of the metric field, which is a scalar 
from the four-dimensional point of view. The trace h 7 ^ acted as an effective source term for the 
graviton propagator in addition to the contribution from the matter source. This complicated 
the tensor structure of the graviton propagator which became momentum dependent. We found 
a solution for the graviton propagator which explicitly revealed its pole structure. We obtained 
infinite towers of massive gravitons and tachyonic ghosts with a discrete mass spectrum. We 
found the contributions from the massive gravitons and tachyonic ghosts in the thin-brane limit 
and showed that the tensor structure and distance dependence of four-dimensional Einstein 
gravity is recovered in this limit. We then analyzed the tensor structure of the momentum 
dependent graviton propagator for a brane of finite thickness. In the small momentum regime, 
the graviton propagator exhibited a .D-dimensional behavior, which was in contrast to the large 
momentum regime (above the critical scale p c (eq. ()44|) ) but well below the inverse brane width A), 
where the contributions from the massive gravitons and tachyonic ghosts conspired to produce 
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a propagator on the brane whose tensor structure and distance dependence was that of four- 
dimensional Einstein gravity. 



[1] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B429 (1998) 263; hep-ph/9803315 



[2] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 (1998) 257; 
|hep-ph? 9804398 

[3] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Rev. D59 (1999) 086004; [hep-ph/9807344 
[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370; |hep-ph /9905221 



[5] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 4690; hep-th/9906064 

[6] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B485 (2000) 208; |hep-t h/0005016, 



[7] G. Dvali, and G. Gabadadze, Phys. Rev. D63 (2001) 065007; hep-th/0008054. 

[8] M. Carena, A. Delgado, J. Lykken, S. Pokorski, M. Quiros and C. E. M. Wagner, Nucl. Phys. B609 

(2001) 499; |hep-ph/0102l72| 
[9] I. I. Kogan, S. Mouslopoulos, A. Papazoglou, G. G. Ross and J. Santiago, Nucl. Phys. B584 (2000) 

313; |hep-ph/99i 2552 



[10] E. Witten, |hep 1 ph/0002297 

[11] R. Gregory V. A. Rubakov and S. M. Sibiryakov, Phys. Rev. Lett. 84 (2000) 5928; 
|hep-th7 0002072 



[12] C. Csaki, J. Erlich and T. J. Hollowood, Phys. Rev. Lett. 84 (2000) 5932; hep-th/0002161. 



[13] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B484 (2000) 112; hep-th/0002190. 



[14] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B484 (2000) 129; hep-th/0003054. 



[15] I. I. Kogan and G. G. Ross, Phys. Lett. B485 (2000) 255; hep-th/0003074 



[16] C. Csaki, J. Erlich, T. J. Hollowood and J. Terning, Phys. Rev. D63 (2001) 065019; 



hep-th/0003076 



[17] I. Giannakis and H. Ren, Phys. Lett. B528 (2002) 133; hep-th/0111127 
[18] H. van Dam and M. Veltman, Nucl. Phys. B22 (1970) 397. 
[19] V. I. Zakharov, JETP Lett. 12 (1970) 312. 

[20] C. Deffayet, G. Dvali, G. Gabadadze and A. I. Vainshtein, Phys. Rev. D65 (2002) 044026; 
|hep-th7 0106001 

14 



[21] A. Lue, Phys. Rev. D66 (2002) 043509; hep-th/0111168 

[22] A. Gruzinov, |astro-pn/0l"l 2246 

[23] M. Porrati, Phys. Lett. B534 (2002) 209; |hep-th/020301^ 

[24] E. Kiritsis, N. Tetradis and T. N. Tomaras, JHEP 0108 (2001) 012; |hep-th/0106050 

[25] G. Dvali, G. Gabadadze, X. Hou and E. Sefusatti, |nep-tn/0i"Tl266 



[26] S. L. Dubovsky, V. A. Rubakov, hep-th/02 12222 



15 



